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EXECUTIVE  SUMMARY 


While  the  general  aviation  industry  increasingly  uses  adhesive  bond  for  airframe  composite 
structures,  two  issues  have  become  particularly  important:  certification  of  adhesives  used  and 
certification  of  adhesive  joints.  A  generally  agreed-upon  design  methodology  for  adhesive- 
bonded  composite  joints,  especially  the  failure  criterion,  is  lacking.  Typically,  three  failure 
modes  of  adhesive-bonded  composite  joints  exist:  (1)  adherend  failure,  (2)  adhesive  failure 
(failure  at  the  adhesive/adherend  interface),  and  (3)  cohesive  failure  (failure  within  the  adhesive 
layer).  A  previous  report  (DOT/FAA/AR-Ol/57,  2001)  included  an  analytical  model 
development  for  adherend  failure  mode.  This  report  describes  the  investigation  of  adhesive  and 
cohesive  failure  modes. 

This  report  includes  two  parts.  The  first  part  comprises  the  derivation  of  an  analytical  model  for 
predicting  adhesive  stress  distribution  within  the  joint  specimens  of  ASTM  D  3165.  Elastic 
orthotropic  adherends  and  elastic-perfectly  plastic  adhesive  were  used  in  the  model  development. 
The  developed  stress  model  was  verified  with  finite  element  models  by  comparing  the  adhesive 
stress  distributions.  Failure  analysis  was  conducted  based  on  three  failure  criteria.  Predicted 
strengths  were  compared  with  test  data.  The  second  part  of  this  report  contains  finite  element 
approaches  for  analyzing  adhesive  and  cohesive  failure  modes.  The  equivalent  plastic  strain  and 
/-integral  were  used  as  the  failure  criteria  for  cohesive  and  adhesive  failure  modes,  respectively. 
Experimental  data  was  used  to  establish  the  values  of  the  equivalent  plastic  strain  and  the  critical 
value  of  /-integral  (/c). 
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1.  INTRODUCTION. 


Advanced  composite  materials  have  been  widely  used  due  to  their  lightweight  and  high- 
corrosion  resistance.  In  many  applications,  bolted  joints  have  been  replaced  by  adhesive-bonded 
joints  because  of  the  weight  penalty  and  corrosion  problems  associated  with  bolted  joints.  Many 
certification-related  issues  become  more  important  as  the  application  of  adhesive-bonded  joints 
gains  in  popularity  in  the  general  aviation  industry.  The  objectives  of  this  study  are  two-fold: 
(1)  to  provide  an  analytical  model  to  predict  the  adhesive  stress  distribution  within  an  adhesive- 
bonded  joint  using  elastic-perfectly  plastic  adhesive  properties  and  to  investigate  the  cohesive 
failure  mode,  and  (2)  to  use  finite  element  methods  to  investigate  the  cohesive  and  adhesive 
failure  modes  of  adhesive-bonded  joints. 

Earlier  studies  on  adhesive-bonded  joints  can  be  found  from  reviewing  papers  by  Kutscha  [1], 
Kutscha  and  Hofer  [2],  Matthews,  et  al.  [3],  and  Vinson  [4].  When  studying  adhesive-bonded 
lap  joints,  the  effects  due  to  the  rotation  of  the  adherends  were  first  taken  into  account  by  Goland 
and  Reissner  [5].  They  introduced  an  equation  to  relate  the  bending  moment  of  the  adherend  at 
the  end  of  the  overlap  to  the  in-plane  loading.  The  basic  approach  of  the  Goland  and  Reissner 
theory  was  based  on  beam  theory  or,  rather,  on  cylindrically  bent-plate  theory,  which  treated  the 
overlap  section  as  a  beam  of  twice  the  thickness  of  the  adherend.  Hart-Smith  [6-9]  published  a 
series  of  papers  regarding  single-lap,  double-lap,  scarf,  and  stepped-lap  joints  involving  a 
continuum  mechanics  model  in  which  the  adherends  were  isotropic  or  anisotropic  elastic,  and  the 
adhesive  was  modeled  as  elastic,  elastic-plastic,  or  bielastic.  Basically,  classical  plate  theory  was 
adopted  during  Hart-Smith’s  derivation.  However,  the  effects  of  transverse  shear  deformation, 
which  has  been  shown  to  be  important  when  the  span-to-depth  ratio  or  the  transverse  shear 
modulus  is  small  (Reissner  [10]  and  Reddy  [11]),  were  not  included  in  either  Goland  and 
Reissner  or  Hart-Smith’s  theories.  Moreover,  edge  effects  were  neglected,  and  adhesive  stresses 
were  assumed  to  be  constant  through  the  thickness  in  most  of  the  analyses  found  in  the  literature. 

Renton  and  Vinson  [12]  used  a  higher  order  formulation  that  includes  the  adherend  transverse 
shear  and  normal  strains  to  analyze  adhesive-bonded  joints.  Their  results  were  compared  with 
photoelastic  experiments  and  show  excellent  agreement,  except  for  local  stress  concentrations  at 
the  edges  of  the  overlap  region.  Yang,  et  al.  [13]  studied  double-lap  composite  joints  under 
cantilevered  bending  and  developed  a  strain  gap  model  to  describe  the  stress-strain  behavior. 
Yang  and  Pang  [14  and  15]  also  derived  analytical  models  for  adhesive-bonded  composite 
single-lap  joints  under  cylindrical  bending  and  tension  based  on  the  laminated  anisotropic  plate 
theory.  Their  approach  correlated  the  asymmetry  of  the  adherend  laminates  as  well  as  the  effects 
due  to  the  transverse  shear  deformation.  Oplinger  [16]  studied  single-lap  joints  with  isotropic 
adherends  and  found  that  many  of  the  predictions  of  the  Goland-Reissner  analysis  are  recovered 
in  the  limit  of  large  adherend-to-adhesive  layer  thickness  ratios,  although  substantial  differences 
from  the  Goland-Reissner  analysis  can  occur  for  relatively  thin  adherends. 

In  1997,  Tong  [17]  studied  the  strength  of  adhesive-bonded  double-lap  composite  joints.  Due  to 
the  fact  that  failure  often  occurs  at  the  resin-fiber  interface  adjacent  to  the  adhesive,  Tong  used  a 
simplified  one-dimensional  model  as  well  as  a  finite  element  model  in  conjunction  with  several 
existing  and  new  interlaminar  failure  criteria  to  predict  the  strength  of  joints.  In  1996,  Adams 
and  Davies  [18]  published  the  results  of  their  nonlinear  finite  element  analysis  on  single-lap 
adhesive-bonded  joints  of  composite/steel  and  composite/aluminum  with  different  taper 


1-1 


arrangements  at  the  edges  of  joints.  An  experimental  study  of  the  effects  of  adherend  layup 
sequence  on  joint  strength  following  ASTM  D  1002  was  published  by  Thomas,  et  al.  [19]  in 
1998. 

Many  researchers  investigated  the  failure  of  adhesive-bonded  joints.  To  name  a  few,  the  linear 
elastic  fracture  mechanics  approach  was  used  by  Chai  [20  and  21]  and  Chai  and  Chiang  [22]. 
They  used  the  Butterfly  and  the  End-Notch  Flexure  adhesive-bond  specimens  to  establish  a  crack 
propagation  criterion  for  adhesive  bonds.  Dynamic  crack  growth  of  bonds  was  investigated  by 
Needleman  and  Rosakis  [23].  The  finite  element  method  has  been  widely  used  to  analyze 
adhesive-bonded  composite  structures  as  well  as  adhesive-bonded  repairs.  Among  these 
researchers,  Charalambides,  et  al.  [24  and  25]  studied  bonded  repairs  under  static  and  fatigue 
loading,  experimentally  as  well  as  using  elastic-plastic  finite  element  models. 

Although  finite  element  analysis  can  solve  many  mechanical  problems  with  different  materials 
and  configurations,  analytical  solutions  are  still  preferred  to  perform  parametric  analyses  such  as 
optimization.  In  the  first  part  of  this  report,  an  analytical  model  is  derived  to  determine  the  stress 
and  strain  distributions  of  adhesive-bonded  composite  single-lap  joints  under  tension.  In  the 
model  development,  the  laminated  anisotropic  plate  theory  is  applied  in  the  derivation  of  the 
governing  equations  of  the  two  bonded  laminates.  The  adhesive  is  assumed  to  be  very  thin,  and 
the  adhesive  stresses  are  assumed  constant  through  the  bondline  thickness.  The  composite 
adherends  are  assumed  linear  elastic  while  the  adhesive  is  assumed  elastic-perfectly  plastic 
following  von  Mises  yield  criterion.  The  entire  coupled  system  is  then  determined  through  the 
kinematics  and  force  equilibrium  of  the  adhesive  and  the  adherends.  The  overall  system  of 
governing  equations  is  solved  analytically  with  appropriate  boundary  conditions.  The  computer 
software.  Maple  V,  is  used  as  the  solution  tool.  Results  from  the  analytical  model  are  verified 
with  finite  element  analysis  using  ABAQUS.  The  three  major  failure  modes  of  adhesive-bonded 
joints  are:  (1)  adherend  failure,  (2)  cohesive  failure,  and  (3)  adhesive  failure.  While  the 
adherend  failure  mode  was  investigated  and  reported  in  the  previous  report  (DOT/FAA/AR- 
01/57,  2001),  the  current  report  includes  cohesive  and  adhesive  failure  modes.  The  analytical 
model  developed  in  this  investigation  is  used  to  predict  the  cohesive  failure  using  experimental 
results  from  ASTM  D  3165-95  “Strength  Properties  of  Adhesives  in  Shear  by  Tension  Loading 
of  Single-Lap-Joint  Laminated  Assemblies”  test  specimen  configurations. 

The  second  part  of  this  report  describes  the  finite  element  analyses  of  test  specimens  specified  in 
both  ASTM  D  3165  and  ASTM  D  5656-95  test  standards.  Both  cohesive  and  adhesive  failure 
modes  are  included  in  the  finite  element  analyses. 
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2.  ANALYTICAL  MODEL  ON  ASTM  D  3165  SPECIMEN. 


As  described  previously,  an  analytical  model  is  developed  to  predict  the  stress  and  strain 
distributions  within  an  adhesive-bonded  joint.  In  the  model,  the  laminated  anisotropic  plate 
theory  describes  the  behavior  of  the  adherends.  The  adhesive  is  assumed  to  be  elastic-perfectly 
plastic  and  to  follow  von  Mises  yield  criterion.  The  entire  system  of  coupled  ordinary 
differential  equations  is  solved  by  using  a  symbolic  solver,  Maple  V  [26].  Finite  element 
analysis  is  conducted  to  confirm  the  adhesive  stress  distributions  obtained  from  the  developed 
model.  Cohesive  failure  is  investigated  using  the  developed  model  and  experimental  results. 

2.1  MODEL  DERIVATION. 

Figure  2-1  shows  the  configuration  of  a  single-lap  joint  following  ASTM  D  3165  configuration. 
The  tensile  load  in  the  x  direction,  defined  as  P,  represents  the  load  per  unit  width.  As  shown 
below,  the  joint  is  divided  into  three  regions  for  convenience  in  the  model  development.  There  is 
an  adhesive  plastic  zone  at  each  edge  of  the  overlap  within  which  the  adhesive  shear  stress 
reaches  its  yield  strength. 
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FIGURE  2-1.  JOINT  CONFIGURATION  AND  COORDINATE  SYSTEMS 
2.1.1  Adherend  Formulation. 

The  generalized  formulas  of  the  adherends  and  adhesive  are  the  same  for  all  three  regions.  The 
upper  and  lower  adherends  are  assumed  to  follow  the  first-order  laminated  plate  theory.  The 
displacement  field  of  the  two  adherends,  u  in  the  x  direction  and  w  in  the  z  direction,  can  be 
written  as  follows: 


(x)  +  zuiffU  (x) 

(2-1) 

(x)  +  ZL  if  O) 

(2-2) 

2-1 


U-  U(Y\ 

w  ~w  W 


(2-3) 


wL  —  wL(x)  (2-4) 

where  superscript  U  and  L  denote  the  upper  and  lower  adherends,  respectively,  superscript  o 
represents  the  mid-plane  displacement,  and  y/  is  the  corresponding  bending  slope.  After 
substituting  equations  2-1  through  2-4  into  the  strain  displacement  relations  together  with  the 
equivalent  modulus  matrices  [A],  [E],  and  [Z>]  for  orthotropic  laminates,  the  normal  stress 
resultant  Nx,  bending  moment  for  unit  width  My,  and  transverse  shear  stress  resultant  Qz  can  be 
obtained  as  follows: 


_  du°(x)  dyr(x) 

Nx  Au  ,  +  Bn 

dx  dx 

(2-5) 

M,  =  B»duy  +  DnJr}X) 

dx  dx 

(2-6) 

_  ,  (  t  .  dw(x)  \ 

a=M55  y 

V  dx  J 

(2-7) 

where  ks  is  the  shear  correction  factor.  In  order  to  establish  the  equations  of  equilibrium,  a  tree- 
body  diagram  of  the  adherends  and  adhesive  is  shown  in  figure  2-2.  The  equations  for  force 
equilibrium  of  the  upper  adherend  are  then  as  follows: 

*3 

i 

II 

^“■8 

(2-8) 

dMUy  „uyu 

dt  ”a  Tr” 

(2-9) 

d(f2_n 

dx  °2 

(2-10) 

where  rxz  and  az  are  the  shear  and  peel  stresses  of  the  adhesive,  respectively,  and  hu  is  the 
thickness  of  the  upper  adherend.  Three  equilibrium  equations  can  be  obtained  for  the  lower 
adherend  in  the  same  fashion. 
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FIGURE  2-2.  FREE-BODY  DIAGRAM  AND  SIGN  CONVENTIONS 
2.1.2  Adhesive  Behavior  in  Elastic  Regions. 

Assuming  a  perfect  bond  between  the  adhesive  and  the  adherend  surfaces,  based  on  the 
kinematics  of  the  adherends,  the  adhesive  strains  are  related  to  the  bottom  surface  of  the  upper 
adherend  and  to  the  top  surface  of  the  lower  adherend.  In  terms  of  the  displacement  field  of  the 
two  adherends,  the  adhesive  strains  can  be  written  as  follows: 


jL  yu 

fh  l  ,  h  — u \ 


7  L  ,U 

2  2 


9^ 

'w'' 

1 _ 1 

<5 

(2-11) 

r) 

(2-12) 

(2-13) 

£z - (wU  -  WL) 

1 


where  rj  is  the  adhesive  thickness.  Under  plane-strain  condition,  the  adhesive  stresses  can  be 
obtained  as  follows: 


g-"o+vXl-*)[(l~y)fo+KglJ 


rju  .  . 

ff'=(l  +  K)(l-2v)fe  +  £j 


ffi=a+7xr>)[l'£-+(1'v)£j 


(2-14) 

(2-15) 

(2-16) 
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and 


Txz  =  GYxz 


(2-17) 


2.1.3  Adhesive  Behavior  in  Plastic  Regions. 

High  adhesive  stresses  at  the  vicinity  of  the  joint  edges  are  expected.  Therefore,  it  is  necessary 
to  assume  a  plastic  zone  near  each  of  the  edges  where  adhesive  yielding  is  considered. 

Based  on  Prandtl-Reuss  material  model,  the  plastic  strain  increment  is  expressed  by  the  flow 
rule: 


id-p-M  (2-18) 

OCTij 

where  sy  and  cr,y  are  the  strain  and  stress  tensors  of  the  adhesive,  respectively,  and  g  is  plastic 
potential.  When  elastic-perfectly  plastic  behavior  following  von  Mises  yield  criterion  is  assumed 
for  the  adhesive,  the  plastic  potential  can  be  assumed  as  follows: 

g-J.-Y,*,  (2-19) 


Therefore,  equation  2-18  can  be  written  as  follows: 

d£%_  _  ^£y_  _  de%_  _  d  7^  _  d  _  d_  7^  _  ^  (2-20) 

Sx  Sy  S  z  Txy  T yz  Tzx 

Also,  for  a  perfectly  plastic  material,  the  equation  can  be  written  as  follows: 

dat  -  2  Ode,  +  KdeaS„  -  gj--f s,  (2-21) 

K 

where  e,y  is  the  deviatoric  strain  tensor,  <5,y  is  the  Kronecker  delta,  k  =  (J2)2,  Ji  is  the  second 
invariant  of  the  deviatoric  stress  tensor,  G  is  the  shear  modulus,  and 


K  = 


E 

3(1 -2v) 


(2-22) 


von  Mises  stress  o'  can  be  expressed  in  terms  of  the  deviatoric  stresses  and  the  shear  stress  as 
follows: 


o'  =  J-(sx + si + si) + 3  t\z 


(2-23) 
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By  assuming  that  the  shear  stress  is  much  greater  than  all  the  three  deviatoric  stress 
components  when  the  adhesive  reaches  the  yield  stress  at  the  edges  of  the  joint,  then 


__  _  O' yield 


(2-24) 


and 


sx~0  sy~0  sz  ~  0  (2-25) 

Based  on  this  assumption  and  equation  2-21,  it  can  be  shown  that  d<Jx  ~  dOy  ~  dcrz  ~  0  after  the 
adhesive  reaches  yield  strength.  Therefore,  the  adhesive  behavior  in  the  plastic  zones  is  as 
follows: 


and 


(l+v)(l-2v) 


[(l-v)£x  +  v£z\ 


(l+v)(l-2v) 


(l+v)(l-2v) 


[v^  +  (l-v)f2] 


(2-26) 

(2-27) 

(2-28) 


_  _  _  O yield 


(2-29) 


This  assumption  is  verified  by  the  finite  element  analysis  using  ABAQUS.  Figure  2-3  shows  the 
adhesive  stress  distribution  in  region  2,  as  specified  in  figure  2-1.  It  can  be  seen  that  the  shear 
stress  is  almost  constant  within  the  adhesive  plastic  zones,  while  the  von  Mises  stress  remains 
constant. 
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FIGURE  2-3.  ADHESIVE  STRESS  DISTRIBUTION  FROM  FEM  (ABAQUS) 

Because  the  adhesive  shear  stress  is  assumed  constant  within  the  plastic  zone,  the  von  Mises 
stress  calculated  based  on  equation  2-23  is  no  longer  a  constant.  A  factor  is  then  applied  to  each 
of  the  stress  components  to  scale  down  the  von  Mises  stress  to  the  value  of  the  yield  stress  ffyieid- 
Due  to  the  symmetry  of  the  joint  system,  the  length  of  the  plastic  zone  in  region  1  {L\p)  is  the 
same  as  the  length  of  the  plastic  zone  in  region  3  (L3p)  and  the  two  plastic  zones  in  region  2  have 
the  same  length  {Lip  =  Lip ’). 

2.1.4  Solution  Procedure. 

The  overall  system  of  governing  equations  includes  42  coupled  second-order  ordinary 
differential  equations  (six  equations  for  each  zone)  with  42  variables  (uoU,  y/u,  wu,  uoL,  ,  and 
WL  for  each  zone).  A  total  of  84  boundary  conditions  are  obtained  at  the  ends  of  each  zone  from 
either  continuity  or  applied  force  conditions,  as  shown  in  figure  2-1.  These  boundary  conditions 
are  listed  in  appendix  A.  The  solution  procedure  starts  with  assumed  plastic  zone  lengths,  L\p 
and  Lip.  Solutions  to  the  differential  equations  are  obtained  by  using  a  symbolic  solver.  Maple  V 
[26],  with  the  Laplace  option.  Once  the  Laplace  option  is  chosen,  Maple  V  performs  forward 
and  inverse  Laplace  transformations  to  obtain  the  exact  solutions  to  the  system  of  ordinary 
differential  equations.  Based  on  the  solutions,  the  von  Mises  (a')  stresses  in  the  plastic  zones  are 
calculated.  The  lengths  of  the  plastic  zones  for  the  next  iteration  are  adjusted  by  multiplying 
(o'/oyieid)2  by  the  plastic  zone  lengths  of  the  previous  iterations.  Iteration  continues  until  the 
calculated  von  Mises  stresses  o'  converge  to  the  value  of  yield  stress  oyieid- 

2.2  RESULTS  AND  DISCUSSION. 

To  demonstrate  the  application  of  the  developed  model,  a  joint  system  as  described  in  figure  2-1 
is  modeled  using  both  the  developed  model  and  finite  element  software  ABAQUS.  In  this 
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illustration,  T300/5208  (Graphite/Epoxy)  with  a  ply  thickness  of  0.25  mm  was  used  for  both 
upper  and  lower  adherends.  Each  adherend  consists  of  12  plies  with  orientation  and  sequence  of 
[03/903]s.  The  engineering  constants  of  T300/5208  are  as  follows: 

E\\  =  181  GPa  E22  =  10.3  GPa  Gi2  =  7.17  GPa  viZ  =  0.28 

For  convenience,  other  mechanical  properties  of  the  adherends  are  assumed  as  £33  =  E22, 
Gn=  G12,  vj3  =  V12,  v23  =  0.35,  and  G23  =  £22/(2+  2v23)  =  3.815  GPa.  Based  on  the  material 
properties  above,  the  mechanical  constants  of  the  adherends  per  unit  width  can  be  given  as 
follows: 


^nU  =  v4iiL  =  374kN/mm  £iiU  =  £nL  =  0 

DuU  =  Dnh  =  394  kNrnrn  A55U  =  A55L  =  18.  2kN/mm 

The  adhesive  is  Metbond  408  with  the  following  properties: 

£  =  0.96  GPa  G  =  0.34  GPa 
v  =  0.41  CTyjeid  =  7.5  MPa 

The  joint  dimensions  include  L\  =  £3  =  80  mm,  L2  =  30  mm,  and  tj  =  0.2  mm.  The  applied 
tensile  load  is  P  =  0.1  kN/mm. 

Figure  2-4  shows  the  von  Mises  stress  distribution  obtained  from  both  the  developed  model  and 
finite  element  model.  A  good  correlation  can  be  seen  from  the  size  of  the  adhesive  plastic  zone 
and  the  stress  level  in  the  entire  overlap  region.  Figures  2-5  and  2-6  show  the  adhesive  peel 
stress  (<rz)  and  shear  stress  distributions,  respectively,  obtained  from  the  developed  model  and  the 
finite  element  model.  It  can  be  seen  that  the  reduced  value  of  shear  stress  at  the  ends  of  the 
overlap  can  be  obtained  following  scaling,  as  previously  described. 

Joint  strength  can  be  predicted  if  a  suitable  joint  failure  criterion  is  established.  Experimental 
results  using  ASTM  D  3165  specimens  with  PTM&W  ES6292  paste  adhesive  from  previous 
investigation  (reported  in  DOT/FAA/AR-Ol/33)  are  used  for  the  establishment  of  a  valid  joint 
failure  criterion.  The  yield  strength  ( oyieid. )  of  PTM&W  ES6292  is  41.37  MPa  based  on  previous 
test  results. 

Several  proposed  criteria  include  the  following:  (1)  hydrostatic  stress  at  certain  locations  such  as 
5%  of  overlap  length  from  the  edges  of  the  joint,  (2)  average  hydrostatic  stress  over  a  portion  of 
the  overlap  length,  and  (3)  average  plastic  strain  over  a  portion  of  the  overlap  length.  The  results 
all  show  good  correlation  with  test  results.  The  index  using  the  hydrostatic  stress  of  the  adhesive 
at  certain  locations  of  the  overlap  can  be  expressed  as  follows: 

Index  1  =  ~(<Tx  +  (7y  +  (Tz)  (2-30) 
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FIGURE  2-4.  ADHESIVE  VON  MISES  STRESS  DISTRIBUTION 


FIGURE  2-5.  ADHESIVE  PEEL  STRESS  DISTRIBUTION 
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x  (mm) 

FIGURE  2-6.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION 
Table  2-1  provides  the  Index i  values  at  a%  of  L2  from  the  overlap  edge  using  test  results. 
TABLE  2-1.  VALUES  OF  INDEXx  FROM  TEST  RESULTS  (MPa) 


A 

0 

2 

3 

4 

5 

Sample  No.  1 

35.65 

28.54 

24.82 

21.10 

17.51 

Sample  No.  2 

36.40 

29.23 

25.44 

21.65 

18.00 

Sample  No.  3 

39.23 

31.44 

27.10 

22.75 

18.48 

Sample  No.  4 

37.16 

29.79 

25.92 

21.93 

18.13 

Sample  No.  5 

39.78 

31.85 

27.44 

22.96 

18.55 

Average 

37.65 

30.17 

26.15 

22.08 

18.13 

Std.  Dev. 

18.91 

1.42 

1.10 

0.77 

0.42 

If  the  average  hydrostatic  stress  within  a  certain  portion  close  to  the  edges  of  the  overlap  is  used 
as  failure  criterion,  Index 2  is  written  as  follows: 

Ilf  ^(Ox  +  ay+Gz)dx 

Index  2  ~  ~ - « -  (2-31) 

ap 

Table  2-2  provides  the  values  of  Index 2  from  test  results. 
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TABLE  2-2.  VALUES  OF  INDEXj  FROM  TEST  RESULTS  (MPa) 


p  =  Overlap  Length 

P  =  Adhesive  Plastic  Region 

A 

5% 

8% 

10% 

30% 

50% 

Sample  No.  1 

31.10 

28.54 

33.92 

30.48 

27.30 

Sample  No.  2 

31.85 

29.23 

34.47 

30.75 

27.30 

Sample  No.  3 

41.23 

31.44 

37.03 

32.54 

28.48 

Sample  No.  4 

32.47 

29.79 

35.16 

31.23 

27.58 

Sample  No.  5 

34.75 

31.85 

37.44 

32.89 

28.61 

Average 

34.28 

30.17 

35.60 

31.58 

27.85 

Std.  Dev. 

17.63 

1.42 

1.56 

1.08 

0.64 

If  the  average  adhesive  plastic  strain  is  used  to  evaluate  failure,  Index3  is  as  follows: 

Index, =  (2-32) 

aP 

Table  2-3  provides  the  values  of  lndex3  from  test  results. 

TABLE  2-3.  VALUES  OF  INDEX3  FROM  TEST  RESULTS  (in/in) 


P  =  Overlap  Length 

P  =  Adhesive  Plastic  Region 

A 

3% 

5% 

8% 

10% 

30% 

50% 

Sample  No.  1 

0.0435 

0.0401 

0.0353 

0.0455 

0.0389 

0.0330 

Sample  No.  2 

0.0480 

0.0444 

0.0395 

0.0496 

0.0423 

0.0358 

Sample  No.  3 

0.0640 

0.0596 

0.0533 

0.0657 

0.0557 

0.0469 

Sample  No.  4 

0.0517 

0.0480 

0.0427 

0.0534 

0.0454 

0.0383 

Sample  No.  5 

0.0675 

0.0628 

0.0562 

0.0692 

0.0585 

0.0492 

Average 

0.0549 

0.0510 

0.0454 

0.0567 

0.0482 

0.0406 

Std.  Dev. 

0.0104 

0.0098 

0.0090 

0.0103 

0.0085 

0.0071 

The  validation  of  the  developed  analytical  model  can  be  justified  based  on  the  reasonable 
correlations  between  the  stress  distributions  obtained  from  the  developed  analytical  model  and 
the  finite  element  model.  However,  due  to  the  fact  that  most  of  the  composite  joints  failed  at  the 
adherends  instead  of  the  adhesive  in  the  experiments  conducted  on  ASTM  D  3165  specimens, 
there  is  not  sufficient  data  to  establish  a  reliable  failure  criterion.  In  tables  2-1  through  2-3, 
results  from  only  five  specimens  are  used  to  establish  the  failure  criterion.  More  data  points 
would  be  necessary  for  the  development  of  a  valid  failure  criterion.  This  would  be  feasible  only 
if  a  ductile  and  relatively  weak  adhesive  is  to  be  used. 
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3.  FINITE  ELEMENT  ANALYSES  OF  ASTM  D  5656  AND  ASTM  D  3165  SPECIMENS. 


3.1  FINITE  ELEMENT  MODELS  OF  ASTM  D  5656  SPECIMENS. 

The  ASTM  D  5656  test  standard  allows  one  to  obtain  the  stress-strain  relationship  of  adhesives. 
From  the  test,  one  can  determine  the  adhesive’s  apparent  shear  modulus,  the  stress  and  strain  at 
the  proportional  limit,  and  the  failure  load  of  the  joint.  The  specimen  is  loaded  at  both  ends, 
while  the  relative  displacement  of  the  two  adherends  is  recorded  using  the  KGR-1  device.  The 
thick  nature  of  the  adherends  causes  them  to  have  a  high  rigidity  in  both  tension  and  bending. 
The  displacement  of  the  adherends  in  the  gauge  section  is  fairly  uniform.  The  ASTM  D  5656 
specimen  geometry  and  dimensions  are  shown  in  figure  3-1. 
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FIGURE  3-1.  ASTM  D  5656  SPECIMEN  GEOMETRY  AND  DIMENSIONS 

To  simulate  the  mechanical  behavior  of  the  joint,  including  stress  distributions  and  failure,  two- 
dimensional  (2-D)  finite  element  models  with  four-node  elements  are  constructed  for  specimens 
with  bondline  thicknesses  of  0.013",  0.04",  0.08",  and  0.12".  The  number  of  elements  of  the 
entire  specimen,  the  number  of  elements  of  the  adhesive,  and  the  element  aspect  ratios  are 
provided  in  table  3-1. 

TABLE  3-1.  NUMBER  OF  ELEMENTS  AND  ASPECT  RATIO  OF  FINITE  ELEMENT 

MODELS  (ASTM  D  5656) 


FE  Model 

Adhesive 

Thickness 

(in) 

Number  of 
Elements  of 
Entire  Model 

Number  of 
Elements  of 
Adhesive 

Number  of 
Elements  Through 
Adhesive 
Thickness 

Largest  Element 
Aspect  Ratio  of 
Adhesive 

D  5656 J 

0.013 

13,420 

400 

8 

5.4 

D  5656  2 

0.040 

10,208 

700 

14 

3.1 

D  5656  3 

0.080 

8,204 

900 

18 

2.1 

D  5656  4 

0.120 

7,864 

880 

22 

2.3 

The  following  techniques  are  used  in  the  finite  element  analyses. 
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3.1.1  Explicit  Direct- Integration  Dynamic  Analysis. 

This  method  is  appropriate  for  quasi-static  analysis  of  large  models  with  nonlinear  material 
properties  and  large  deformation.  The  explicit  direct-integration  dynamic  analysis  is  available  in 
AB  AQU S/Explicit.  It  is  based  on  the  implementation  of  explicit  central  difference  integration 
together  with  the  use  of  diagonal  element  mass  matrices.  The  explicit  schemes  obtain  values  for 
dynamic  quantities  at  t  +  At  directly  from  the  available  values  at  time  t,  which  is  efficient  for 
each  time  increment.  Therefore,  a  large  number  of  small  time  increments  are  performed 
effectively.  However,  the  central  difference  operator  is  conditionally  stable.  The  stability  limit 
(the  largest  time  step)  is  approximately  equal  to  the  time  for  an  elastic  wave  to  move  across  the 
smallest  element  in  the  model.  In  the  present  analysis,  the  time  increments  are  determined 
automatically  by  ABAQUS/Explicit  within  the  stability  limit.  Moreover,  the  mass  scaling 
method  is  used  in  the  calculation  to  expedite  the  simulation  process. 

3.1.2  Nonlinear  Analysis  for  Both  Adhesive  and  Adherends. 

Both  geometry  and  material  nonlinearity  are  included  in  the  finite  element  simulation.  To 
include  material  nonlinearity,  the  von  Mises  yield  criterion  with  isotropic  hardening  and  the 
associated  plastic  flow  are  used,  which  can  be  expressed  in  equations  3-1  and  3-2. 

f  =  4T1-k(e„)  =  0  (3-1) 

ds'=dk^~  (3-2) 

d<T„ 

where  J2  is  the  second  invariant  of  the  deviatoric  stress  tensor;  k  is  the  value  defining  the 
occurrence  of  material  yielding,  which  is  a  function  of  effective  plastic  strain  £p  for  isotropic 

hardening  materials;  dA  is  a  positive  scalar  factor  of  proportionality;  /  is  von  Mises  yield 
function;  and  cr  andff  are  stress  tensor  and  plastic  strain  tensor,  respectively.  The  effective 

plastic  strain  £p  is  defined  as  follows: 


The  incremental  stress-strain  relation  can  be  generally  expressed  as  follows: 


d(7,j  —  y,  (CiikI  +Cgti  )d£kI 

A:  ,7=1, 2, 3 


(3-4) 


where  i,  j  =1,  2,  3,  the  incremental  stress  tensor  da^  and  strain  tensors  dsy  can  be  written  in 
vector  forms  as  follows: 


{daf  =  {d<Jxd<7yd  (Tzd  T^d  Zad  } 


(3-5) 
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{de}T  =  {de^s^e^e^dsjej 


(3-6) 


The  tensor  of  elastic  modulus  Qju  and  C?jkl  are  expressed  in  the  following  matrix  forms: 
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where,  sx,  sy,  sz,  Sxy,  s^,  and  are  the  reduced  stress  deviatoric  tensors,  and  G  and  K  are  the 
shear  and  bulk  moduli,  respectively.  Hp  and  the  effective  stress  ae  are  expressed  as  follows: 


H 


p 


(3-9) 


(3-10) 


Actually,  Hp  is  the  slope  of  the  uniaxial  true  stress/true  plastic  strain  curve,  which  is  required  by 
ABAQUS  as  input  data. 

The  material  used  for  the  adherends  is  2024-T351  aluminum.  Its  material  properties  are 
simulated  using  an  elastic-perfectly  plastic  model  with  Young’s  modulus  of  10.6  Msi,  Poisson’s 
ratio  of  0.3,  and  yield  stress  of  47  ksi. 

PTM&W  ES6269  paste  adhesive  is  used  in  the  finite  element  analyses  using  the  elastic-plastic 
material  model  with  isotropic  hardening.  In  the  elastic  range,  the  Young’s  modulus  of  the 
adhesive  is  275  ksi  and  the  Poisson’s  ratio  is  0.3.  The  postyielding  behavior  required  by 
ABAQUS  is  a  plastic  strain/true  stress  curve.  Because  the  material  data  is  usually  available  in 
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values  of  nominal  stress  and  strain,  it  is  necessary  to  convert  the  nominal  stress/strain  relation  to 
the  true  stress/strain  relation  before  it  can  be  used  by  ABAQUS.  The  true  strain  or  logarithmic 
strain  e and  true  stress  crare  defined  as  follows: 


£- 


=  ln 


(3-11) 


F 


(3-12) 


where  /  is  length,  l0  is  the  original  length,  F  is  the  applied  force,  and  A  is  the  cross-sectional  area. 
According  to  the  definitions,  the  relations  among  true  stress,  true  strain,  nominal  stress,  and 
nominal  strain  are  as  follows: 


a +0  (3‘13) 

£=M1+£U)  (3-14) 

The  plastic  strain  is  obtained  by  subtracting  the  elastic  strain,  which  is  defined  as  the  value  of 
true  stress  divided  by  the  Young’s  modulus,  from  the  total  strain,  or  the  following  equation: 

£pl  =£-crl  E  (3-15) 

To  convert  material  test  data  to  the  required  ABAQUS  input,  equations  3-13  and  3-14  are  used  to 
convert  the  nominal  stress  and  nominal  strain  to  true  stress  and  true  strain.  Plastic  strain  is 
determined  from  true  stress  and  true  strain  using  equation  3-15.  The  converted  stress  and  strains 
are  shown  in  table  3-2.  Although  only  a  slight  difference  between  the  nominal  and  true  values 
exists  at  small  strains,  there  is  a  significant  difference  at  larger  strains.  Therefore,  it  is  important 
to  provide  the  proper  stress-strain  data  to  ABAQUS  if  large  strains  occur  in  the  simulation. 
Figure  3-2  shows  the  nominal  stress  versus  nominal  strain,  true  stress  versus  true  strain,  and  true 
stress  versus  plastic  strain  curves. 

In  fact,  the  adhesive  material  data  in  table  3-2  is  obtained  by  trial  calculations  using  finite 
element  models  and  compared  with  experimental  results  using  ASTM  D  5656  standard  test  with 
a  bondline  thickness  of  0.013".  As  shown  in  figure  3-3,  the  experimental  results  correlate  well 
with  the  finite  element  simulation  using  the  material  data  in  table  3-2. 
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TABLE  3-2.  POSTYIELDING  BEHAVIORS  OF  PTM&W  ES6269  ADHESIVE 


Nominal  Stress 

Gnom 

Nominal  Strain 

£nom 

True  Stress 

G 

True  Strain 

s 

Plastic  Strain 
il 

0.0000 

0.0000 

0.0000 

0.0000 

3.2601 

0.0122 

3.3000 

0.0121 

0.0000 

3.9356 

0.0163 

4.0000 

0.0162 

0.0015 

4.4172 

0.0187 

4.5000 

0.0185 

0.0020 

5.0053 

0.0231 

5.1210 

0.0228 

0.0040 

5.4399 

0.0294 

5.6000 

0.0290 

0.0084 

5.7416 

0.0380 

5.9600 

0.0373 

0.0154 

5.9111 

0.0456 

6.1810 

0.0446 

0.0219 

6.0265 

0.0548 

6.3570 

0.0533 

0.0300 

6.1121 

0.0686 

6.5320 

0.0664 

0.0424 

6.1251 

0.0799 

6.6150 

0.0769 

0.0526 

6.0954 

0.1049 

6.7350 

0.0997 

0.0750 

5.9845 

0.1330 

6.7810 

0.1249 

0.1000 

5.7212 

0.1913 

6.8160 

0.1750 

0.1500 

5.4352 

0.2523 

6.8070 

0.2250 

0.2000 

Strain 

FIGURE  3-2.  MATERIAL  PROPERTIES  OF  PTM&W  ES6292  PASTE  ADHESIVE 
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FIGURE  3-3.  COMPARISON  BETWEEN  EXPERIMENTAL  RESULTS  (ASTM  D  5656) 
AND  FE  SIMULATION  USING  MATERIAL  PROPERTIES  DATA  IN  TABLE  3-2 

3.1.3  Displacement  Control  at  the  End  of  Adherend. 

To  simulate  the  actual  loading  process  in  the  experiment,  displacement  boundary  conditions  are 
used  at  the  loading  end  of  the  specimen.  The  displacement  increases  linearly  along  the  loading 
direction  at  a  rate  of  0.05  in/min,  which  is  the  same  as  in  the  experiment. 

3.1.4  Half  Model. 

According  to  the  loading  condition  and  geometry  of  the  specimen,  only  half  of  the  specimen  is 
modeled  with  appropriate  boundary  constraints.  As  figure  3-4  shows,  point  O  is  the  geometry 
center  of  the  specimen  and  line  A-A  is  perpendicular  to  the  loading  direction.  Along  line  A-A, 
the  displacements  of  two  symmetric  points  B  and  B'  with  respect  to  O  must  satisfy  the  following 
constraints: 

uBx+u°=  0,  uBy+uB  =  0  (3-16) 

where  uBx  and  uB  are  the  displacements  along  the  x  direction  of  points  B  and  B'  and  uy  and 

uy  are  the  displacements  along  y  direction  of  points  B  and  B',  respectively.  The  constraint 

equations  are  established  for  all  the  nodes  on  line  A-A  in  the  finite  element  model.  By 
implementing  these  constraint  equations  into  the  finite  element  model  together  with  a  zero- 
displacement  boundary  condition  at  point  O,  only  half  of  the  specimen  is  modeled  to  reduce  the 
computation  time.  The  right  half  of  the  specimen  is  selected  in  the  present  study. 


V 
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FIGURE  3-4.  SYMMETRIC  BOUNDARY  CONDITIONS  AT  CENTER  OF 

ASTM  D  5656  SPECIMEN 


3.1.5  Shear  Failure  Model. 

According  to  the  experiment  results,  failure  of  ASTM  D  5656  specimen  occurs  when  the 
adhesive  reaches  large  plastic  strains.  Therefore,  failure  models  for  ductile  materials  are 
appropriate  to  simulate  ASTM  D  5656  specimen  behavior.  ABAQUS/Explicit  features  a  shear 
failure  criterion,  which  is  used  in  the  present  study  to  simulate  the  failure  of  adhesive  (cohesive 
failure  mode).  The  shear  failure  model  is  based  on  the  value  of  the  equivalent  plastic  strain  at 
element  integration  points.  Failure  is  assumed  to  occur  when  the  damage  parameter  exceeds  one. 
The  damage  parameter  0) is  defined  as  follows: 


\  ) 


where  ep>  is  the  equivalent  plastic  strain  defined  as 


V  3  ij= 1,2,3 


(3-17) 


(3-18) 


and  is  designated  as  PEEQ  in  ABAQUS/Explicit.  The  £p  is  plastic  strain  rate  tensor,  As pl  is 
the  increment  of  PEEQ,  and  ef  is  the  plastic  strain  at  failure.  The  summation  in  equation  3-17 

is  performed  over  all  increments  in  the  analysis.  When  the  shear  failure  criterion  is  reached  at  an 
element  integration  point,  all  the  stress  components  are  set  to  zero  and  the  material  fails  at  that 
point.  If  the  material  fails  at  all  integration  points  of  an  element,  the  element  is  removed  from 
the  mesh.  In  the  present  study,  the  first-order  reduced-integration  2-D  solid  elements  are  used. 
The  value  of  ef  is  chosen  as  0.2  according  to  trial  calculations  and  experiment  results. 
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3.2  FINITE  ELEMENT  MODELS  OF  ASTM  D  3165  SPECIMENS. 


Two-dimensional  finite  element  models  with  eight-node  elements  are  used  to  simulate  ASTM  D 
3165  specimens  with  bondline  thicknesses  of  0.013",  0.04",  0.08",  and  0.12",  as  shown  in  figpre 
3-5.  The  number  of  elements  of  the  entire  specimen,  the  number  of  elements  of  the  adhesive, 
and  the  element  aspect  ratios  are  given  in  table  3-3. 


FIGURE  3-5.  ASTM  D  3165  SPECIMEN  GEOMETRY  AND  DIMENSIONS 

TABLE  3-3.  NUMBER  OF  ELEMENTS  AND  ASPECT  RATIO  OF  FINITE  ELEMENT 

MODELS  (ASTM  D  3165) 


4 


FE  Model 

Adhesive 

Thickness 

(in) 

Number  of 
Elements  of 
Entire  Model 

Number  of 
Elements  of 
Adhesive 

Number  of 
Elements  Through 
Adhesive 
Thickness 

Largest  Element 
Aspect  Ratio  of 
Adhesive 

D  31651 

0.013 

6,020 

580 

8 

5.1 

D  31652 

0.040 

4,284 

860 

14 

3.5 

D  3 1 65 _ 3 

0.080 

3,872 

920 

18 

2.8 

D  31654 

0.120 

3,944 

1,120 

22 

2.5 

Only  half  of  the  ASTM  D  3165  specimen  is  modeled  due  to  the  symmetry  of  both  loading  and 
specimen  geometry  (as  presented  previously  for  the  finite  element  models  of  ASTM  D  5656 
specimens).  Displacement  control  is  applied  at  the  loading  end  of  model.  The  material 
properties  of  both  the  adherend  and  adhesive  are  the  same  as  those  used  for  ASTM  D  5656 
specimens.  The  geometric  and  material  nonlinearities  are  also  included  and  the  nonlinear 
equilibrium  equations  are  solved  with  Newton’s  method  in  ABAQU S/Standard.  The  automatic 
increment  scheme  is  used  because  the  increment  is  based  upon  computational  efficiency. 

Static  analysis  is  used  for  two  reasons:  (1)  the  nonlinearity  of  the  material  and  the  deformation 
in  the  adhesive  of  ASTM  D  3165  specimens  are  not  as  severe  as  those  of  ASTM  D  5656  at 
failure  and  (2)  methods  of  fracture  mechanics,  such  as  the  7-integral  criterion,  critical  crack 
opening  displacement  criterion,  can  be  applied  with  static  analysis. 

The  techniques  of  fracture  mechanics  analysis  associated  with  the  finite  element  model  are 
described  as  follows. 
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3.2.1  /-Integral. 


The  /-integral  is  applied  in  quasi-static  fracture  analysis  in  ABAQUS/Standard  to  characterize 
the  energy  release  rate  associated  with  crack  growth.  Considering  an  arbitrary  counter-clockwise 
path  (/)  around  the  tip  of  a  crack,  as  illustrated  in  figure  3-6,  the  /-integral  is  as  follows: 


(3-19) 


where  w  is  the  strain  energy  density;  uj,U2,  and  us  are  the  displacement  vector  components;  ds  is 
the  length  increment  along  the  contour  F,  and  Tj,  T2,  and  T3  are  components  of  the  traction 
vector.  The  traction  is  a  stress  vector  normal  to  the  contour.  In  other  words,  /;,  T2,  and  T3  are 
the  normal  stresses  acting  at  the  boundary  if  a  free-body  diagram  on  the  material  inside  of  the 
contour  is  constructed.  Theoretically,  the  values  of  /-integral  should  be  independent  of  the  path 
selected.  However,  the  calculated  /-integral  based  on  different  rings  may  vary  because  of  the 
approximate  nature  of  finite  element  solutions.  Therefore,  at  least  two  contours  are  used  for  /- 
integral  calculation  according  to  the  ABAQUS/Standard  User’s  Manual.  Based  on  the 
convergence  test,  the  /-integral  on  the  fourth  ring  of  elements  around  the  crack  front  shows 
appropriate  stability. 


FIGURE  3-6.  ARBITRARY  CONTOUR  AROUND  TIP  OF  A  CRACK 
3.2.2  Modeling  the  Crack  Tip  Singularity. 

Because  the  failure  mode  of  the  aluminum-adherend  ASTM  D  3 165  specimen  is  considered  to  be 
the  fracture  between  the  adherend  and  adhesive,  i.e.,  the  adhesive  failure  mode,  a  precrack  must 
be  built  into  the  finite  element  model.  The  /-integral  calculated  on  the  fourth  ring  is  used  as  the 
index  forjudging  the  onset  of  the  crack.  To  improve  the  accuracy  of  the  /-integral  and  the  stress 
and  strain  calculations,  the  singularity  at  the  crack  tip  is  modeled  as  suggested  in 
ABAQUS/Standard  Manual.  Eight-node  isoparametric  elements  must  be  used  and  one  side  of 
the  element  is  collapsed  so  that  all  three  nodes  a,  b,  and  c  have  the  same  geometric  location  (on 
the  crack  tip)  to  fulfill  the  mesh  singularity  requirement,  as  depicted  in  figure  3-7.  During  this 
collapsing  process,  the  middle  node  remains  in  the  middle  and  nodes  a,  b,  and  c  are  allowed  to 
move  independently.  In  this  case,  the  1/r  (r  is  the  distance  from  the  crack  tip)  singularity  in 
strain  is  created,  which  is  accurate  enough  for  /-integral  calculation  in  the  present  study 
according  to  trial  calculations.  The  typical  mesh  with  the  precrack  is  shown  in  figure  3-8.  The 
length  of  the  precrack  is  approximately  3%  of  the  length  of  the  overlap  for  all  models. 
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FIGURE  3-7.  COLLAPSED  TWO-DIMENSIONAL  ELEMENT 


FIGURE  3-8.  TYPICAL  FOCUSED  MESH  FOR  FRACTURE  MECHANICS  EVALUATION 
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3.3  RESULTS  AND  DISCUSSION  OF  ASTM  D  5656  SPECIMENS. 


It  is  difficult  to  gain  more  detailed  information  from  the  actual  ASTM  D  5656  experiment  other 
than  the  apparent  shear  stress-strain  curves  and  apparent  shear  strength.  The  failure  mode  is 
estimated  by  examining  the  failed  specimen,  which  does  not  reflect  what  actually  occurs  during 
the  test.  However,  by  using  the  finite  element  method,  the  mechanical  behavior  of  ASTM  D 
5656  specimens  during  the  loading  process  can  be  simulated.  The  results  from  finite  element 
simulation  can  provide  more  detailed  information  regarding  deformation,  stress  and  strain 
distributions,  damage  process,  and  damage  mode  of  the  specimen.  This  information  is  important 
for  understanding  the  damage  mechanism  of  such  joints  and  in  realizing  the  load  carrying 
capability,  which  provides  practical  guidance  in  joint  design. 

3.3.1  Deformation  Analysis  (ASTM  D  56561. 

Detailed  deformation  analysis  in  the  overlap  area  of  the  specimen  is  very  important.  It  is  the  first 
step  towards  understanding  the  mechanical  behavior  of  the  joint  during  loading.  For  ASTM  D 
5656  specimens  with  different  bondline  thickness,  figures  3-9  through  3-12  show  the  calculated 
displacements,  without  applying  any  failure  criterion,  along  the  top  and  bottom 
adhesive/adherend  interfaces  under  each  respective  experimental  failure  load.  In  these  figures,  ux 
and  uy  are  the  displacements  along  the  load  direction  and  the  direction  perpendicular  to  the  load 
direction,  respectively,  as  shown  in  figure  3-4.  It  can  be  seen  that  the  displacements  are  linear 
with  respect  to  x,  for  all  four  bondline  thicknesses.  The  angles  of  rotation  of  the  adherends  in  the 
overlap  region  for  different  bondline  thicknesses  under  the  respective  failure  load  are  listed  in 
table  3-4.  The  small  angles  indicate  that  bending  does  not  contribute  much  to  the  adherend 
deformation  within  the  overlap  due  to  the  large  adherend  thickness. 

TABLE  3-4.  ANGLE  OF  ROTATION  OF  OVERLAP  AREA  OF 
ASTM  D  5656  SPECIMENS  AT  FAILURE  LOAD 


Bondline 
Thickness,  t, 

(in) 

Experimental  Apparent 
Adhesive  Shear  Strain 
at  Failure 

Experimental 
Failure  Load 
(kip) 

Rotation  Angle  of 
Overlap  at  Failure 
(FE  Model) 

0.013 

0.30 

1.41 

0.152E 

0.040 

0.20 

1.34 

0.1 82E 

0.080 

0.12 

1.26 

0.208E 

0.120 

0.065 

1.05 

0.212E 

Figures  3-13  and  3-14  show  the  through-thickness  average  adhesive  shear  and  peel  strain 
distributions  along  the  overlap  for  specimens  with  different  bondline  thicknesses  when  each 
respective  experimental  failure  load  is  applied  to  the  finite  element  model.  The  through¬ 
thickness  average  adhesive  shear  strain  in  figure  3-13  is  calculated  by  dividing  the  difference  in 
displacement  in  the  x  direction  between  the  upper  and  lower  adhesive/adherend  interfaces  by  the 
bondline  thickness.  Similarly,  the  through-thickness  average  adhesive  peel  strain  is  calculated 
by  dividing  the  difference  of  displacements  in  the  transverse  direction  (y  direction).  It  can  be 
seen  from  figure  3-13  that  the  through-thickness  average  shear  strain  is  fairly  uniform  within  the 
overlap  and  that  larger  average  adhesive  shear  strains  exist,  at  failure,  in  joints  with  thinner 
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bondlines.  Figure  3-14  indicates  that  the  through-thickness  average  adhesive  peel  strain  does 
vary  along  the  overlap  length  for  specimens  with  a  thinner  bondline,  especially  at  the  end  of  the 
overlap.  However,  the  through-thickness  average  peel  strain  is  small  compared  to  the  through¬ 
thickness  average  shear  strain.  It  should  be  noted  that  these  through-thickness  average  adhesive 
strains  cannot  fully  determine  the  actual  stress  distributions  because  both  the  ffee-surface 
boundary  conditions  and  stress  concentration  at  the  edges  of  the  overlap  have  great  effects  on  the 
adhesive  stress  distributions. 


FIGURE  3-9.  DISPLACEMENT  OF  ASTM  D  5656  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.013"  AT  FAILURE  LOAD  (1.41  kips) 
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FIGURE  3-10.  DISPLACEMENT  OF  ASTM  D  5656  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.040"  AT  FAILURE  LOAD  (1.34  kips) 
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FIGURE  3-11.  DISPLACEMENT  OF  ASTM  D  5656  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.080"  AT  FAILURE  LOAD  (1.26  kips) 
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FIGURE  3-12.  DISPLACEMENT  OF  ASTM  D  5656  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.120"  AT  FAILURE  LOAD  (1.05  kips) 
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FIGURE  3-13.  THROUGH-THICKNESS  AVERAGE  ADHESIVE  SHEAR  STRAIN  WITHIN 
OVERLAP  AT  EACH  RESPECTIVE  FAILURE  LOAD  (ASTM  D  5656) 


FIGURE  3-14.  THROUGH-THICKNESS  AVERAGE  ADHESIVE  PEEL  STRAIN  WITHIN 
OVERLAP  AT  EACH  RESPECTIVE  FAILURE  LOAD  (ASTM  D  5656) 

3.3.2  Stress  Analysis  (ASTM  D  5656). 

The  adhesive  stress  distributions  are  quite  complex  in  the  ASTM  D  5656  specimen.  Generally 
speaking,  shear  is  the  main  mechanism  which  transfers  load  from  one  adherend  to  the  other. 
However,  due  to  the  load  eccentricity  of  the  adherends,  both  the  adhesive  shear  stress  and  peel 
stress  exist  at  the  ends  of  the  overlap  and  play  an  important  role  in  failure  assessment. 
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Theoretically,  the  adhesive  stresses  in  the  ASTM  D  5656  specimen  can  be  fully  determined  by 
the  displacement  boundary  conditions  at  the  upper  and  lower  adhesive/adherend  interfaces  and 
the  two  free  adhesive  surfaces  at  the  ends  of  the  overlap.  Because  the  tensile  rigidity  (EA  value) 
of  the  adherends  is  relatively  high,  the  displacement  along  the  load  direction  at  the  two 
adhesive/adherend  interfaces  is  rather  uniform.  Therefore,  the  main  factor  which  creates  the 
complexity  of  the  adhesive  stresses  is  the  free  adhesive  surfaces  at  the  ends  of  the  overlap. 

Figure  3-15  shows  five  paths  that  are  along  the  load  direction  and  evenly  distributed  through  the 
bondline  thickness.  The  adhesive  shear  stress  distribution  under  a  tensile  load  of  800  lb  for  the 
five  paths  is  shown  in  figures  3-16  through  3-19.  It  can  be  seen  that  the  adhesive  shear  strain 
decreases  to  zero  at  the  free  adhesive  surface  (Paths  2,  3,  and  4  at  x  =  0.1875")-  The  portion  of 
the  overlap  where  the  adhesive  shear  stress  deviates  from  the  near-constant  value  is  obviously 
larger  when  the  adhesive  free  surface  at  the  ends  of  the  overlap  is  larger,  i.e.,  in  the  case  of 
thicker  bondlines.  The  adhesive  free  surface  actually  causes  this  nonuniformity  of  the  shear 
stress  distribution.  Figure  3-20  shows  the  five-path  average  adhesive  shear  stress  distribution 
based  on  the  four-bondline  thicknesses.  Similarly,  thicker  bondlines,  or  larger  adhesive  free 
surfaces  at  the  ends  of  the  overlap,  affect  a  larger  portion  of  the  overlap.  As  a  result,  under  the 
same  applied  load,  in  order  to  compensate  for  the  larger  affected  overlap  area  by  the  free  surface, 
the  adhesive  shear  stress  in  the  middle  portion  is  higher  for  specimens  with  larger  bondline 
thicknesses.  In  the  actual  ASTM  D  5656  tests,  the  apparent  adhesive  shear  stress  is  calculated  by 
dividing  the  applied  load  by  the  overlap  area.  The  larger  adhesive  shear  stress  at  the  middle 
overlap  portion  in  the  finite  element  model  is  reflected  as  larger  adhesive  shear  strain  in  the 
actual  ASTM  D  5656  tests.  Therefore,  the  higher  adhesive  shear  strain  in  the  middle  portion  of 
the  overlap  for  the  thicker  bondline  explains  the  fact  that  thicker  bondlines  result  in  lower 
adhesive  shear  moduli  when  the  ASTM  D  5656  test  is  conducted.  Table  3-5  shows  the  finite 
element  model  simulated  ASTM  D  5656  test  results.  It  is  obvious  that  lower  adhesive  shear 
moduli  are  calculated  for  larger  bondline  thicknesses. 

It  should  be  noted  that  if  the  bondline  becomes  very  thin  (thinner  than  0.01"  for  adhesives  with 
Gi  100  ksi,  [27])  the  shear  modulus  obtained  from  ASTM  D  5656  tests  is  actually  larger  than  the 
true  value.  This  is  because  the  effect  of  nonuniformity  of  the  adherend  deformation  is  no  longer 
negligible  when  the  adhesive  is  very  thin.  Actually,  the  shear  stress  distribution  in  figure  3-20 
has  already  shown  this  trend  where  the  five-path  average  adhesive  shear  stress  distribution  starts 
peaking  up  at  the  edge  of  the  overlap  for  the  case  of  0.013"  bondline  thickness. 


FIGURE  3-15.  DEFINITION  OF  FIVE  PATHS  IN  ADHESIVE  LAYER 


3-15 


Shear  Stress  (ksi) 


FIGURE  3-16.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.013"  BONDLINE 

THICKNESS,  ASTM  D  5656) 


FIGURE  3-17.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.04"  BONDLINE 

THICKNESS,  ASTM  D  5656) 
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FIGURE  3-18.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.08"  BONDLINE 

THICKNESS,  ASTM  D  5656) 


FIGURE  3-19.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.12"  BONDLINE 

THICKNESS,  ASTM  D  5656) 
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FIGURE  3-20.  FIVE-PATH  AVERAGE  ADHESIVE  SHEAR  STRESS  DISTRIBUTION 

(ASTM  D  5656) 

TABLE  3-5.  ADHESIVE  SHEAR  MODULUS 


Bondline  Thickness  (in) 

0.013 

0.04 

0.08 

0.12 

Apparent  Shear  Modulus,  Gastm  (ksi), 
Finite  Element  Simulated 

99.28 

94.73 

89.99 

86.02 

Deviation* 

-6.1% 

-10.4% 

-14.9% 

-18.7% 

_  ...  G  ASTM  G input 

Note:  Deviation  - - 


Where  GASTM  is  the  calculated  adhesive  shear  modulus  using  finite  element  models  to  simulate  the  actual  ASTM  D 
5656  test  and  G,„pu(=105.77  ksi  is  the  material  property  used  in  the  finite  element  model. 

Figure  3-21  shows  the  development  of  shear  stress  distribution  along  path  3  (as  defined  in  figure 
3-15),  while  the  displacement  at  the  end  of  the  joint  gradually  increases  (as  shown  in  figure 
3-22).  When  the  load  nearly  reaches  the  failure  load  of  the  joint,  the  adhesive  has  already 
entered  its  plastic  state.  After  the  adhesive  plastic  strain  has  initiated,  both  the  load  and  the 
adhesive  stress  rise  more  slowly  even  though  the  joint-end  displacement  remains  at  the  same 
speed. 
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FIGURE  3-21.  DEVELOPMENT  OF  SHEAR  STRESS  DISTRIBUTION  ALONG  PATH  3 
(0.12"  BONDLINE  THICKNESS,  ASTM  D  5656) 


FIGURE  3-22.  LOAD  VERSUS  JOINT-END  DISPLACEMENT  (0.12"  BONDLINE 

THICKNESS,  ASTM  D  5656) 

3.3.3  Failure  Mode  and  Failure  Criterion  (ASTM  D  56561. 

As  stated  in  the  ASTM  D  5656  test  standard,  the  purpose  of  this  test  is  to  determine  the  stress- 
strain  relationship  of  the  adhesive.  However,  understanding  the  failure  mode  and  failure 
criterion  can  assist  in  design  and  analysis  of  such  joints. 
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The  three  major  failure  modes  of  adhesive-bonded  joints  are  as  follows:  (1)  cohesive  failure, 
(2)  adhesive  failure,  and  (3)  adherend  failure.  Because  ASTM  D  5656  specimens  use  aluminum 
adherends,  their  failure  tends  to  be  cohesive  and/or  adhesive  failures.  Cohesive  failure  is 
characterized  by  the  failure  of  the  adhesive  itself.  In  the  case  of  cohesive  failure,  the  adhesive  is 
left  on  both  bonded  surfaces  of  the  adherends.  However,  sometimes  it  is  difficult  to  judge 
cohesive  failure  for  a  specimen  because  of  the  random  propagation  of  the  failure  surface  inside 
the  adhesive.  The  failure  surface  might  develop  inside  the  adhesive,  propagate  into  the 
adhesive/adherend  interface,  and  continue  to  develop  along  the  interface.  It  is  also  worthwhile  to 
note  that  even  though  cohesive  failure  is  the  failure  mode  of  the  specimens,  the  apparent  strength 
of  the  joint,  which  is  obtained  by  dividing  the  failure  load  by  the  adhesive  area  (overlap  area), 
cannot  be  regarded  as  merely  material  dependent.  Their  values  are  closely  related  to  the 
configuration  of  the  joint  as  well. 

Figure  3-23  shows  ASTM  D  5656  test  results  using  four  different  bondline  thicknesses.  The 
load/apparent  adhesive  shear  strain  curves  for  different  bondline  thicknesses  (t)  all  have  a  ramp 
section  and  a  descending  or  near-horizontal  section.  In  the  ramp  section,  the  apparent  shear 
modulus  deviates  from  one  bondline  thickness  to  another,  as  described  previously.  After 
yielding,  the  apparent  shear  strain  at  failure  is  several  times  larger  for  the  case  of  a  thin  bondline 
than  for  the  case  of  a  thick  bondline.  This  means  that  when  the  joint  fails,  there  would  be  a  much 
larger  plastic  strain  accumulated  inside  the  adhesive  layer  of  thinner  bondlines.  The  observation 
that  all  the  curves  are  fairly  smooth,  even  at  the  turning  sections,  implies  that  the  damage 
developing  process  inside  the  adhesive  is  stable.  No  staggered  impetuous  damage  occurs  during 
the  loading  process.  Therefore,  it  is  reasonable  to  assume  that  the  failure  is  caused  by  the 
accumulation  of  plastic  strain  in  the  adhesive. 

Based  on  this  assumption,  the  Shear  Failure  Criterion  featured  in  ABAQUS/Explicit,  as 
described  in  “3.1  (5)  Shear  Failure  Model,”  is  applied  to  investigate  the  failure  load  of  ASTM  D 
5656  specimens.  To  further  demonstrate  the  dominating  factors  of  joint  failure,  peel  stress,  shear 
stress,  and  PEEQ  (equivalent  plastic  strain)  at  two  critical  locations  within  the  adhesive  layer  are 
analyzed  for  specimens  with  different  bondline  thicknesses.  Figure  3-24  shows  the  critical 
locations  where  stress  concentration  occurs  during  the  loading  process.  Location  1  is  at  the 
comer  of  the  adhesive  layer,  which  is  the  first  adhesive  element  along  the  lower  bondline  at  the 
free-boundary  surface.  Location  2  is  where  the  actual  failure  is  initiated  in  the  finite  element 
model  with  a  bondline  thickness  of  0.013".  It  is  located  at  the  upper  interface  between  the 
adhesive  and  adherend  and  is  about  15%  of  the  overlap  length  away  from  the  free  boundary  of 
the  adhesive. 


3-20 


FIGURE  3-23.  EXPERIMENTAL  LOAD  VERSUS  APPARENT  ADHESIVE  SHEAR 

STRAIN  CURVES  (ASTM  D  5656) 


FIGURE  3-24.  CRITICAL  LOCATIONS  INSIDE  ADHESIVE  LAYER 

Figures  3-25  to  3-27  show  the  shear  stress,  peel  stress,  and  PEEQ  development  at  the  two  critical 
locations,  respectively.  Figures  3-25  and  3-26  show  that  it  is  impossible  to  find  some  critical 
stresses  as  an  index  of  failure,  because  the  adhesive  materials  at  the  critical  locations  are  in  large 
plastic  states  in  which  the  stress  has  almost  no  change.  The  stress  criteria  cannot  reflect  the 
actual  damage  state  in  this  case.  Different  from  stresses,  PEEQ  always  increases  during  loading 
until  the  critical  value  of  0.2,  as  shown  in  figure  3-27.  The  critical  value  of  0.2  is  carefully 
chosen  by  trial  computation  in  order  to  correlate  the  failure  loads  and  the  corresponding  apparent 
shear  strains  at  failure  of  all  specimens  in  the  experiments.  Note  that  the  initial  damage  location 
of  the  specimen  with  bondline  thickness  of  0.013"  is  different  from  other  specimens.  As  shown 
in  figure  3-27,  its  PEEQ  increases  faster  at  location  1  than  location  2  initially,  but  the  PEEQ  at 
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Shear  Stress 


location  2  grows  much  faster  and  finally  exceeds  the  PEEQ  at  location  1  before  reaching  the 
critical  value.  Finally,  the  initial  damage  occurs  at  location  2. 


FIGURE  3-25.  SHEAR  STRESS  AT  CRITICAL  LOCATIONS 
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FIGURE  3-26.  PEEL  STRESS  AT  CRITICAL  LOCATIONS 


FIGURE  3-27.  PEEQ  AT  CRITICAL  LOCATIONS 


Figure  3-28  shows  the  calculated  apparent  shear  stress-strain  curves  with  the  equivalent  plastic 
shear  strain  as  failure  criterion.  The  curves  for  thin  bondlines  correlate  well  with  the 
experimental  results  shown  in  figure  3-29.  The  main  deviation  is  that  the  failure  loads  and  the 
corresponding  apparent  shear  strains  at  failure  of  thick  bondline  (0.08"  and  0.12")  from  the  finite 
element  models  are  higher  than  those  from  the  experiments.  The  reason  might  be  that  the 
localized  effects  due  to  the  free  boundary  are  more  severe  in  the  adhesive  of  the  specimens  with 
thicker  bondlines.  The  PEEQ  at  the  critical  locations  accumulates  quickly  in  the  loading  process 
while  damage  initiates  in  the  elements  at  the  critical  locations  after  PEEQ  reaches  the  critical 
value.  However,  damage  does  not  develop  as  quickly  among  the  elements  near  the  damaged 
areas  compared  to  what  actually  occurred  in  experiments.  This  is  because  damage  initiation  and 
damage  development  in  the  finite  element  models  should  be  judged  with  different  criteria  or,  at 
least,  the  same  criterion  with  different  parameters.  Usually,  the  criterion  for  damage  propagation 
is  lower  than  the  criterion  for  damage  initiation.  In  the  present  study,  only  the  shear  failure 
criterion  with  a  constant  parameter  is  used  forjudging  the  entire  damage  process,  which  includes 
damage  initiation  and  development.  Damage  in  the  specimens  with  thin  bondlines  show  that 
plastic  strain  develops  more  evenly  in  the  adhesive.  Once  the  elements  at  the  critical  location 
become  critical,  the  neighboring  elements  of  the  damaged  elements  reach  the  critical  value 
almost  simultaneously.  Therefore,  the  calculation  results  seem  more  accurate  for  the  specimens 
with  thin  bondlines. 
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FIGURE  3-28.  CALCULATED  APPARENT  SHEAR  STRESS-STRAIN  CURVES  FOR 

ASTM  D  5656  SPECIMENS 


FIGURE  3-29.  EXPERIMENTAL  APPARENT  SHEAR  STRESS-STRAIN  CURVES  FOR 

ASTM  D  5656  SPECIMENS 
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Figures  3-30  to  3-33  show  the  process  of  PEEQ  development  inside  the  adhesive  with  four 
different  bondline  thicknesses,  0.013",  0.04",  0.08",  and  0.12"  respectively.  Each  picture 
describes  the  PEEQ  distribution  at  four  stages.  The  first  stage  represents  the  first  appearance  of 
the  plastic  state  inside  the  adhesive  layer  (the  element  with  little  lighter  color  at  the  lower  right 
comer).  The  second  stage,  shown  in  figures  3-31  through  3-33,  is  the  moment  of  the  first 
appearance  of  damaged  elements,  i.e.,  PEEQ  reaches  0.2  for  this  element  at  that  time.  Because 
there  is  no  element  damage  occurring  before  the  final  fracture  for  the  specimen  with  a  bondline 
thickness  of  0.013",  the  second  stage,  in  figure  3-30,  gives  the  distribution  of  PEEQ  at  the 
moment  that  the  entire  adhesive  reaches  plastic  state.  The  other  two  stages  represent  the 
moments  right  before  and  after  the  final  fracture,  respectively. 

From  these  figures,  it  can  be  seen  that  the  element  at  location  1  is  the  first  to  reach  a  plastic  state 
during  the  loading  process,  indicating  that  location  1  has  the  most  serious  stress  concentration. 
Thereafter,  plastic  deformation  gradually  develops  from  location  1  and  also  from  the  upper  right 
comer  in  the  adhesive.  For  the  specimens  with  bondline  thicknesses  of  0.12",  0.08",  and  0.04", 
PEEQ  develops  more  quickly  at  location  1  and  results  in  final  fracture  inside  the  adhesive.  For 
the  specimen  with  the  0.013"  thick  bondline,  PEEQ  develops  quickly  at  location  1  at  the 
beginning  and  then  slows  down.  Meanwhile,  PEEQ  at  location  2  begins  to  increase  much  faster 
and  finally  leads  to  the  initiation  of  damage.  These  pictures  also  show  that  the  damage 
developments  after  the  initiation  of  damage  take  place  more  rapidly  in  the  adhesive  with  thin 
bondlines  than  in  those  with  thick  bondlines.  As  shown  in  figure  3-30,  the  final  fracture  for  the 
specimen  with  the  0.013"  bondline  thickness  occurs  almost  instantly  after  PEEQ  at  location  2 
reaches  the  critical  value.  For  other  specimens  with  thicker  bondlines,  the  damage  gradually 
develops  from  the  initiation  of  damage  to  the  final  fracture  of  the  whole  model,  which  can  be 
seen  from  the  development  of  apparent  shear  strain  between  stages  2  and  3  in  figures  3-31 
through  3-33. 

As  shown  in  figure  3-29,  the  failure  of  joints  with  thick  bondlines  (0.082"  and  0.120")  occurs 
before  the  development  of  large  plastic  deformation.  Therefore,  the  ./-integral  might  be  a  more 
suitable  failure  criterion  than  the  PEEQ  criterion  for  the  cases  of  thick  bondlines.  To  obtain 
more  accurate  predictions  of  failure  loads  for  specimens  with  thick  bondlines,  the  ./-integral 
failure  criterion  is  used  as  an  alternative.  The  experimental  results  are  shown  in  figure  3-29,  and 
the  load  versus  ./-integral  curves  from  the  finite  element  model  are  shown  in  figure  3-34.  Table 
3-6  lists  the  experimental  and  predicted  joint  failure  loads  where  the  critical  J-integral  value  (, Jc ) 
is  taken  as  4  psi-in.  It  can  be  seen  that  the  /-integral  criterion  is  also  a  good  prediction  for  the 
failure  of  the  ASTM  D  5656  joint  with  thicker  bondlines. 
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(3).  Apparent  Shear  Strain  0*2107 
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<4)  Fracture 

Note:  The  scale  coordinates  factor  in  through^ thickness 
direction. is  4,  along  loading  direction  is  1*. 


FIGURE  3-30.  PEEQ  DISTRIBUTION  IN  ADHESIVE  WITH  BONDLINE  THICKNESS  OF 

0.013"  AT  FOUR  TYPICAL  STAGES 
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(2)  Apparent  Shear  Strain  0*1436 
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<3>  Apparent;  Shear  Strain  0*3136 
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<4)  Fracture 


Note t The  scale  coordinates  factor  in  through-thickness 
direction  is  1*5,  along  loading  direction  is  1* 


FIGURE  3-31.  PEEQ  DISTRIBUTION  IN  ADHESIVE  WITH  BONDLINE  THICKNESS  OF 

0.040"  AT  FOUR  TYPICAL  STAGES 
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< 21)  Apparent  Shear  Strain  0*1380 
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<4>  Fracture 

Note:  The  scale  coordinates  factor  In  through-th i ckness 
direction  is  1^  along  loading  direction  is  lt 


FIGURE  3-32.  PEEQ  DISTRIBUTION  IN  ADHESIVE  WITH  BONDLINE  THICKNESS  OF 

0.080"  AT  FOUR  TYPICAL  STAGES 
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Note:  The  scale  coordinates  factor  in  through-thickness 
direction  is  t,  along  loading  direction  is  1* 


FIGURE  3-33.  PEEQ  DISTRIBUTION  IN  ADHI 

0.120"  AT  FOURTY] 
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FIGURE  3-34.  ./-INTEGRAL  FOR  ASTM  D  5656  OF  SPECIMENS  WITH  FOUR 

BONDLINE  THICKNESSES 

TABLE  3-6.  EXPERIMENTAL  AND  PREDICTED  FAILURE  LOADS  BY  /-INTEGRAL 

CRITERION  (ASTM  D  5656) 


Adhesive 
Thickness  (in) 

Experimental 
Failure  Load  (lb) 

Predicted  Failure 
Load  (lb) 

Deviation  (%) 

0.013 

1,405 

1,413 

0.57 

0.040 

1,338 

1,294 

-3.2 

0.080 

1,256 

1,116 

-11 

0.120 

1,049 

987 

-6.0 

3.4  RESULTS  AND  DISCUSSION  OF  ASTM  D  3165  SPECIMENS. 

The  aim  of  the  ASTM  D  3165  standard  test  is  to  obtain  comparative  shear  strengths  of  adhesive 
based  on  the  specimen  configuration  specified  in  the  standard.  The  simulation  results  in  this 
section  will  be  helpful  in  further  understanding  the  failure  mechanism  and  to  determine  the 
dominating  factors  contributing  to  joint  strength. 

3.4.1  Deformation  Analysis  I  ASTM  D  3 1 65). 

Similarly  to  the  analysis  of  ASTM  D  5656  specimens,  the  displacements  along  the  transverse 
direction  (y  direction,  as  shown  in  figure  3-4)  and  along  the  load  direction  at  the  top  and  bottom 
adhesive/adherend  interfaces,  the  through-thickness  average  peel  strain  and  shear  strain  and  the 
rotation  angles  of  the  overlap  areas  are  analyzed  to  obtain  an  overall  picture  of  how  the 
specimens’  configuration  affect  displacement  fields. 
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Figures  3-35  to  3-38  show  the  displacements  of  the  top  and  bottom  adhesive/adherend  interfaces 
in  the  load  direction  and  the  direction  perpendicular  to  the  load  direction  for  specimens  with 
different  bondline  thicknesses  loaded  with  each  corresponding  failure  load.  Compared  with 
those  of  ASTM  D  5656  specimens  in  figures  3-9  to  3-12,  note  that  the  specimens  in  both  the  two 
standard  tests  have  clearly  different  displacement  characteristics.  The  displacements  of  ASTM 
D  3 165  specimens  in  the  transverse  direction  are  more  than  ten  times  larger  than  those  of  ASTM 
D  5656  specimens  with  all  four  bondline  thicknesses.  This  is  because  the  bending  stiffness  of 
the  adherends  of  ASTM  D  3165  specimens  is  much  lower  than  those  of  ASTM  D  5656 
v  specimens.  However,  displacements  along  the  load  direction  at  failure  are  not  as  large  as  those 

of  ASTM  D  5656  specimens,  because  the  ASTM  D  3165  specimens  have  different  failure 
mechanisms  and  their  failure  occurs  before  the  appearance  of  large  plastic  deformation  in  the 
adhesive  layer.  In  contrast,  there  are  large  plastic  deformations  in  the  adhesive  of  ASTM  D  5656 
specimens  at  failure. 


FIGURE  3-35.  DISPLACEMENT  OF  ASTM  D  3165  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.013"  AT  FAILURE  LOAD  (1.40  kips) 
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Displacement  (in) 


FIGURE  3-36.  DISPLACEMENT  OF  ASTM  D  3165  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.040"  AT  FAILURE  LOAD  (1.10  kips) 


FIGURE  3-37.  DISPLACEMENT  OF  ASTM  D  3165  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.080"  AT  FAILURE  LOAD  (0.677  kips) 
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FIGURE  3-38.  DISPLACEMENT  OF  ASTM  D  3165  SPECIMEN  WITH  BONDLINE 
THICKNESS  OF  0.120"  AT  FAILURE  LOAD  (0.510  kips) 

Because  displacements  are  generally  linear  functions  of  x,  the  angle  of  rotation  of  the  overlap 
can  be  estimated.  As  shown  in  table  3-7,  the  angle  of  rotation  for  ASTM  D  3165  specimens  is 
about  ten  times  larger  than  those  of  ASTM  D  5656.  Also,  because  the  adherends  at  the  notches 
have  plastic  deformation,  permanent  rotation  exists  after  the  load  is  removed. 

TABLE  3-7.  ANGLE  OF  ROTATION  OF  OVERLAP  AREA  OF  ASTM  D  3165  SPECIMENS 

AT  FAILURE  LOAD 


Bondline 

Thickness 

(in) 

Failure  Load 
(kips) 

Rotation  Angle  of 
Bondline 

0.013 

1.40 

2.022° 

0.040 

1.10 

2.437° 

0.080 

0.677 

2.043° 

0.120 

0.510 

2.153° 

The  angles  of  rotation  are  large  enough  to  change  the  stress  distribution  inside  the  adhesive  in  the 
deformed  shape.  However,  because  the  deformation  due  to  bending  of  the  adherends  is  mainly 
located  at  the  notches,  the  stresses  at  the  vicinity  of  the  edges  of  the  overlap  are  influenced 
significantly. 

Figure  3-39  shows  the  through-thickness  adhesive  shear  strain  distribution  for  ASTM  D  3165 
specimens,  without  applying  any  failure  criterion,  for  different  bondline  thicknesses  at  each 
respective  experimental  failure  load.  It  can  be  seen  that  the  through-thickness  adhesive  shear 
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strains  for  bondline  thicknesses  of  0.04"  and  0.013"  vary  greatly  along  the  load  direction. 
Similarly,  the  through-thickness  peel  strains  shown  in  figure  3-40  change  greatly  along  the  load 
direction,  especially  for  the  specimen  with  0.013"  bondline  thickness.  Compared  with  the 
through-thickness  average  adhesive  shear  and  peel  strain  of  ASTM  D  5656  specimens,  as  given 
in  figures  3-13  and  3-14,  it  is  obvious  that  the  ASTM  D  5656  specimens  have  approximately 
twice  the  shear  strain  as  ASTM  D  3165  specimens,  but  ASTM  D  3165  specimens  have  much 
greater  adhesive  peel  strain.  These  results  suggest  that  both  the  peel  and  shear  stress  play 
important  rolls  in  the  failure  assessment  of  ASTM  D  3 165  specimens. 


FIGURE  3-39.  THROUGH-THICKNESS  AVERAGE  ADHESIVE  SHEAR  STRAIN  WITHIN 
OVERLAP  AT  EACH  RESPECTIVE  FAILURE  LOAD  (ASTM  D  3165) 


FIGURE  3-40.  THROUGH-THICKNESS  AVERAGE  ADHESIVE  PEEL  STRAIN  WITHIN 
OVERLAP  AT  EACH  RESPECTIVE  FAILURE  LOAD  (ASTM  D  3165) 
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3.4.2  Stress  Analysis  (ASTM  D  3165). 


As  mentioned  in  the  previous  stress  analysis  of  ASTM  D  5656  specimens,  the  adhesive  shear 
stress  is  the  major  mechanism  that  transfers  the  load  from  one  adherend  to  the  other.  The 
adhesive  shear  stress  distribution  determines  the  mechanical  behavior  of  the  joint  and  is  closely 
related  to  joint  failure.  Figures  3-41  to  3-44  show  the  shear  stress  distributions  of  specimens 
with  four  different  bondline  thicknesses  under  a  load  of  800  lbs.  The  five  curves  in  each  figure 
represent  the  five  paths  in  the  adhesive,  as  defined  in  figure  3-15.  The  five-path  average  shear 
stress  distribution  for  each  bondline  thickness  is  shown  in  figure  3-45. 

It  can  be  seen  that  shear  stress  distributions  of  ASTM  D  3165  specimens  vary  dramatically 
compared  with  those  of  ASTM  D  5656  specimens.  The  shear  stress  at  the  middle  of  the  overlap 
is  low,  while  the  shear  stress  near  the  end  of  the  overlap  is  high.  Because  the  adherends  of 
ASTM  D  3165  specimens  are  thinner  than  ASTM  D  5656  specimens,  their  stiffness  is  much 
lower,  which  leads  to  the  nonuniformity  of  shear  stress  distribution.  Stress  nonuniformity  results 
in  lower  joint  strength  because  damage  occurs  in  the  high  stress  area,  while  the  material  at  low 
stress  area  has  not  been  fully  used  for  its  load-carrying  capability. 

Another  important  feature  is  the  large  ffee-boundary  effect  of  ASTM  D  3165  specimens. 
Because  of  the  thinner  adherends  of  the  ASTM  D  3165  specimens  compared  with  ASTM  D  5656 
specimens,  ASTM  D  3165  specimens’  adhesive  has  larger  free  surfaces  in  relation  to  adherend 
thickness,  where  the  adhesive  stress  must  equal  zero. 


FIGURE  3-41.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.013"  BONDLINE 

THICKNESS,  ASTM  D  3 1 65) 
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- Path  4 

-  -  -Path  5 


FIGURE  3-42.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.040"  BONDLINE 

THICKNESS,  ASTM  D  3165) 


- Path  4 

-  -  -  -Path  5 


FIGURE  3-43.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.080"  BONDLINE 

THICKNESS,  ASTM  D  3165) 


FIGURE  3-44.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  (0.120"  BONDLINE 

THICKNESS,  ASTM  D  3165) 


FIGURE  3-45.  FIVE-PATH  AVERAGE  ADHESIVE  SHEAR  STRESS  DISTRIBUTION 

(ASTM  D  3165) 
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Figure  3-46  indicates  the  development  of  shear  stress  during  the  loading  process.  Figure  3-47 
shows  the  development  of  force  while  the  displacement  at  the  loading  end  increases.  Both 
figures  show  the  results  obtained  from  the  specimen  with  bondline  thickness  of  0.12".  It  can  be 
seen  that  the  trend  of  shear  stress  distribution  at  different  loads  are  quite  similar.  By  comparing 
figure  3-47  with  figure  3-22,  it  can  be  seen  that  the  end  joint  displacement  of  ASTM  D  3165 
specimens  is  much  greater  than  those  in  ASTM  D  5656  specimens  at  the  same  load  because 
ASTM  D  3165  specimens  have  much  thinner  adherends. 
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FIGURE  3-46.  DEVELOPMENT  OF  SHEAR  STRESS  DISTRIBUTION  ALONG  PATH  3 
(0.12"  BONDLINE  THICKNESS,  ASTM  D  3165) 


FIGURE  3-47.  LOAD  VERSUS  JOINT-END  DISPLACEMENT  (0.12"  BONDLINE 

THICKNESS,  ASTM  D  3165) 
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3.4.3  Failure  Mode  and  Failure  Criterion  (ASTM  D  3165). 


Based  on  the  experimental  failure  loads  of  the  specimen  (table  3-7)  and  the  load-displacement 
curves  shown  in  figure  3-48,  ASTM  D  3165  specimens  fail  before  the  adhesive  yield  stress  has 
fully  developed  in  the  overlap.  In  other  words,  the  failure  mechanism  of  ASTM  D  3165 
specimens  appears  to  be  different  from  the  failure  mechanism  of  ASTM  D  5656  specimens.  The 
failure  of  ASTM  D  5656  specimens  occurs  after  the  accumulation  of  a  large  plastic  strain. 
However,  ASTM  D  3165  specimens  fail  before  the  load-displacement  curve  levels  off.  This 
suggests  that  the  failure  of  ASTM  D  3165  specimens  might  begin  at  some  critical  location, 
propagate  quickly,  and  finally  rupture.  The  brittle  failure  of  ASTM  D  3165  specimens  also 
suggests  the  possibility  of  using  a  fracture  mechanics  approach,  such  as  the  /-integral,  to  analyze 
their  failure. 


FIGURE  3-48.  LOAD  VERSUS  JOINT-END  DISPLACEMENT  (ASTM  D  3165) 

According  to  the  finite  element  models,  the  critical  location  of  stress  is  at  the  lower  right  comer 
of  the  adhesive,  as  shown  in  figure  3-49.  Figures  3-50  to  3-53  show  the  curves  of  adhesive  shear 
stress,  adhesive  peel  stress,  PEEQ,  and  /-integral  versus  applied  load,  respectively,  at  the  critical 
location.  The  respective  values  of  adhesive  shear  stress,  peel  stress,  PEEQ,  and  /-integral  under 
experimental  failure  loads  are  listed  in  table  3-8.  Clearly,  the  peel  stress,  shear  stress,  and  PEEQ 
cannot  be  used  as  the  failure  index  of  ASTM  D  3 165  specimens.  Based  on  trial  calculation,  if  an 
appropriate  value  of  /-integral  is  chosen  as  the  failure  criterion  of  models  with  a  precrack,  the 
failure  loads  calculated  are  close  to  those  from  experimental  data.  As  listed  in  table  3-8,  a  fairly 
good  prediction  of  failure  loads  can  be  achieved  when  the  critical  /-integral  value  (/c)  is  chosen 
as  4  psi-in.  As  shown  in  table  3-9,  once  Jc  is  chosen  as  4  psi-in,  the  predicted  failure  loads  are 
compared  with  the  experimental  results  for  the  ASTM  D  3165  specimen  with  four  different 
bondline  thicknesses. 
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Locat  ion- 1 


FIGURE  3-49.  CRITICAL  LOCATION  INSIDE  ADHESIVE  LAYER  (ASTM  D  3165) 
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FIGURE  3-50.  ADHESIVE  SHEAR  STRESS  VERSUS  LOAD  AT  LOCATION  1 

(ASTM  D  3165) 


3-' 


FIGURE  3-52.  PEEQ  VERSUS  LOAD  AT  LOCATION  1  (ASTM  D  3165) 
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FIGURE  3-53.  /-INTEGRAL  VERSUS  LOAD  AT  LOCATION  1  (ASTM  D  3165) 
TABLE  3-8.  CORRESPONDING  VALUES  AT  FAILURE  LOADS  FOR  ASTM  D  3165 


Adhesive 

Thickness 

(in) 

Experimental 
Failure  Load 
(lb) 

Shear 

Stress 

(ksi) 

Peel 

Stress 

(ksi) 

PEEQ 

/-integral 

(psi-in) 

0.013 

1,405 

3.88 

8.37 

0.29 

5.1 

0.040 

1,099 

3.64 

8.87 

0.088 

4 

0.080 

677 

3.36 

8.63 

0.05 

4 

0.120 

510 

1.81 

6.91 

0.095 

3.44 

TABLE  3-9.  COMPARISON  OF  PREDICTED  FAILURE  LOADS  WITH  EXPERIMENTAL 
RESULTS  FOR  ASTM  D  3165  USING  Jc  =  4  psi-in 


Adhesive 

Thickness 

(in) 

Experimental 
Failure  Load 
(lb) 

Predicted  Failure 
Load 
(lb) 

Deviation 

(%) 

0.013 

1,405 

1,269 

-9.7 

0.040 

1,099 

934 

-15 

0.080 

677 

678 

0.15 

0.120 

510 

552 

8.2 
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3.5  PARAMETERS  AFFECTING  JOINT  PERFORMANCE. 


3.5.1  Adherend  Thickness. 

With  the  same  approach  used  to  simulate  ASTM  D  3165  specimens,  finite  element  models  with 
double  and  triple  the  original  adherend  thickness  of  ASTM  D  3165  specimens  (0.064")  are  used 
to  show  the  effects  of  adherend  thickness  on  the  joint  strength.  The  calculated  /-integral  values 
as  functions  of  applied  load  for  specimens  with  double  and  triple  adherend  thicknesses  are  shown 
in  figures  3-54  and  3-55,  respectively.  When  using  4  psi-in  as  the  critical  /-integral  value  (/c), 
the  calculated  failure  load  increases  as  the  adherend  thickness  increases,  as  shown  in  figure  3-56. 
Furthermore,  it  can  be  seen  from  figure  3-56  that  the  failure  load  decreases  when  the  adhesive 
becomes  thicker. 


0  0.5  1  1.5  2 

Load  (kip) 


FIGURE  3-54.  /-INTEGRAL  FOR  ASTM  D  3165  CONFIGURATION  WITH  DOUBLE 

ADHEREND  THICKNESS 
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3.5.2  Adherend  Stiffness. 


The  calculated  7-integral  value  versus  applied  load  for  ASTM  D  3165  specimens  with  different 
adherend  moduli  is  shown  in  figure  3-57.  If  4  psi-in  is  used  as  the  critical  value  for  /-integral  as 
the  joint  failure  criterion,  the  relationship  between  the  failure  load  and  adherend  Young’s 
modulus  is  shown  in  figure  3-58.  It  can  be  seen  that  with  the  same  joint  geometry,  a  joint  with 
steel  adherends  would  be  stronger  than  one  with  aluminum  adherends.  With  the  same  approach 
and  the  same  joint  failure  criterion,  figures  3-59  and  3-60  show  the  effects  of  adherend  modulus 
on  joint  strength  of  ASTM  D  5656  specimens.  Because  the  ASTM  D  5656  specimens  already 
have  thick  adherends  with  large  bending  and  tensile  rigidities,  increasing  the  adherend  modulus 
does  not  contribute  to  higher  joint  strength. 


0  0.5  1  1.5  2 

Load  (kip) 


FIGURE  3-57.  /-INTEGRAL  VERSUS  LOAD  FOR  WITH  DIFFERENT  ADHEREND 
YOUNG’S  MODULI  (0.013"  BONDLINE  THICKNESS,  ASTM  D  3165) 
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FIGURE  3-58.  FAILURE  LOAD  VERSUS  ADHEREND  YOUNG’S  MODULUS  (0.013" 

BONDLINE  THICKNESS,  ASTM  D  3165) 


0  0.5  1  1.5  2 

Load  (kip) 

FIGURE  3-59.  /-INTEGRAL  VERSUS  LOAD  FOR  WITH  DIFFERENT  ADHEREND 
YOUNG’S  MODULI  (0.013"  BONDLINE  THICKNESS,  ASTM  D  5656) 
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FIGURE  3-60.  FAILURE  LOAD  VERSUS  ADHEREND  YOUNG’S  MODULUS 
(0.013"  BONDLINE  THICKNESS,  ASTM  D  5656) 

3.5.3  Adherend  Plasticity. 

In  the  finite  element  simulation  conducted  in  this  research,  elastic-perfectly  plastic  behavior  is 
assumed  for  the  aluminum  adherends  of  both  ASTM  D  5656  and  D  3165  specimens.  This  is 
reasonable  for  ASTM  D  5656  specimens  because  the  large  thickness  of  ASTM  D  5656 
aluminum  adherends  results  in  minimal  yielding  in  the  adherend.  In  order  to  validate  the  use  of 
elastic-perfectly  plastic  material  properties  of  aluminum  adherend  for  ASTM  D  3165  specimens, 
comparisons  are  made  between  finite  element  models  with  plastic  and  elastic-perfectly  plastic 
aluminum  properties.  Figures  3-61  and  3-62  show  the  calculated  /-integral  versus  applied  load 
with  two  material  models  for  ASTM  D  3165  joints  with  0.04"  and  0.08"  bondline  thicknesses, 
respectively.  The  fact  that  almost  no  distinction  exists  between  the  two  curves  shown  in  figures 
3-61  and  3-62,  especially  at  the  4  psi-in  7-integral  level,  confirms  the  appropriateness  of  using 
elastic-perfectly  plastic  material  model  for  aluminum  adherends. 
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FIGURE  3-61.  COMPARISON  BETWEEN  DIFFERENT  MATERIAL  MODELS  FOR 
ALUMINUM  ADHERENDS  (0.04"  BONDLINE  THICKNESS,  ASTM  D  3165) 
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FIGURE  3-62.  COMPARISON  BETWEEN  DIFFERENT  MATERIAL  MODELS  FOR 
ALUMINUM  ADHERENDS  (0.08"  BONDLINE  THICKNESS,  ASTM  D  3165) 
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4.  CONCLUSION. 


An  analytical  model  for  adhesive-bonded  single-lap  joints  under  tension  was  developed  based  on 
the  first-order  laminated  plate  theory  for  the  adherend  and  on  elastic-perfectly  plastic  behavior 
for  the  adhesive.  Based  on  the  developed  model,  closed-form  solutions  for  stress  distributions  of 
both  the  laminates  and  the  adhesive  were  obtained.  The  proposed  closed-form  solutions  were 
found  to  correlate  with  the  solutions  obtained  through  finite  element  models.  Several  joint 
failure  criteria  were  proposed  and  evaluated  with  experimental  data. 

Extensive  finite  element  modeling  was  performed  to  analyze  both  ASTM  D  5656  and  D  3165 
specimen  behavior.  With  the  displacement  and  stress  fields  obtained  through  finite  element 
calculation,  phenomena  from  experimental  observations,  such  as  different  shear  moduli  obtained 
from  ASTM  D  5656  specimens  with  different  bondline  thickness,  were  explained.  The  proposed 
failure  criteria,  the  equivalent  adhesive  plastic  strain  and  ./-integral,  were  used  effectively  to 
predict  the  joint  strength  of  ASTM  D  5656  and  D  3165  specimens.  The  calculated  joint 
strengths  correlate  well  with  the  experimental  results.  The  effects  on  joint  strength  due  to 
bondline  thickness,  adherend  thickness,  and  adherend  modulus  were  also  investigated. 
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APPENDIX  A— BOUNDARY  CONDITIONS 
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